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Abstract 

Planar L-loop maximally helicity violating amplitudes in Af = 4 supersymmetric 
Yang-Mills theory are believed to possess the remarkable property of satisfying iteration 
relations in L. We propose a simple new method for studying the iteration relations for 
four-particle amplitudes which involves the use of certain linear differential operators and 
eliminates the need to fully evaluate any loop integrals. We carry out this procedure in 
explicit detail for the two-loop amplitude and argue that this method can be used to prove 
the iteration relations to all loops up to polynomials in logarithms. 



January 2006 



Contents 

[l] Introduction |l| 

|2] Advanced Introduction and Motivation [3| 

pT| The Two-Loop Calculation [5] 

|3X| The ±(mP) 2 Term | 

12] The Mj 2) Term | 



|331 Fixing The Ambiguity |T| 



|4~1 Higher Loops 14 

[4. 1 .| General Constructions 15 

[4.2.| Specific Conjectures 16 

|5~] Summary and Outlook 17 



1. Introduction 

Maximally supersymmetric Af = 4 Yang-Mills (SYM) theory possesses remarkably rich 
mathematical structure which has been the subject of intense investigation over the past 
several years. One motivation for much of this work is the AdS/CFT correspondence |]J, 
which asserts that the strongly-coupled SYM theory admits an equivalent formulation as 
gravity in AdSs, thereby opening a new window for studying quantum gravity. A com- 
plementary motivation, which has seen a dramatic resurgence following (see for a 
review) , is the desire to explore the mathematical structure of Yang-Mills perturbation the- 
ory and to exploit that structure to aid the calculation of scattering amplitudes. Scattering 
amplitudes are the basic building blocks which enter into the calculation of experimentally 
measured processes. As the LHC comes on line in the next couple of years there will be 
increasing pressure to bring the theoretical uncertainties in QCD calculations, especially 
at higher loops, under control. 

Optimistically, we hope that the rich structure of Yang-Mills perturbation theory and 
the simplicity of the strongly-coupled theory expected from AdS/CFT are two sides of the 
same coin, and that we might in some cases be able to see some hint of the structure which 
enables perturbation theory to be resummed to match onto AdS/CFT. 

One intriguing step in this direction has been the study of iterative relations amongst 
planar maximally helicity violating (MHV) loop amplitudes in dimensionally regulated 
(d = 4 — 2e) J\f = 4 SYM In [|J Anastasiou, Bern, Dixon and Kosower suggested 

that two-loop MHV amplitudes obey the iteration 

AfJ 2) (e) = \ (M«(e)) 2 + f^\e)M^{2e) - ^ + 0(c), (1.1) 
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where we use M« (e) = A^n\e)/An^ to denote the ratio of the L-loop n-particle color- 
stripped partial amplitude to the corresponding tree-level amplitude, and 

/(2) (e) = ± m _ e) _ m) = _ (C(2)+C(3)e + C(4)e 2 + ... )5 ^ {x) = ±\nT{x). (1.2) 

We refer to ( |1 . 1| ) for general n as the ABDK conjecture. In section 4 we review an ex- 
ponentiated ansatz conjectured in [|5| which generalizes the iterative structure ( p. . 1|) for 
MHV amplitudes to all loops. These conjectures were motivated by studies of the infrared 
|],[7|,|8| and collinear [f3||9|4T2"|| behavior of multiloop amplitudes, where similar basic iterative 
structures appear. 

The problem of calculating L > 1 amplitudes and then verifying these iterative struc- 
tures involves numerous extremely difficult technical challenges. So far the only iteration 
relations which have been explicitly verified are the two-loop ABDK relation ( |1 . 1|) for n = 4 
particles (the amplitude M\ (e) was originally calculated in |L3[ using unitarity based 



methods [ |I4| - PU| ), an d its three-loop generalization [[|, again just for n = 4, where the 
calculation of (e) by Bern, Dixon and Smirnov through O(e ) is an impressive tour 
de force. Even at two loops the n > 4 amplitudes are not yet known. For some progress 
in this direction see JTB| , |2"T|] . 

In this paper we introduce new methods for studying these iteration relations for n = 4 
particle amplitudes. We provide a direct check of the form of the ABDK relation at two 
loops and then show that these methods can easily be generalized to higher L, where they 
should greatly facilitate the problem of verifying conjectured iteration relations. We begin 
by expressing each amplitude appearing in ( |1 . 1|) in a familiar integral form. We then apply 
a linear differential operator U^' which is chosen in such a way that the terms in ( |1 . 1| ) can 
be nicely combined into a single integral which is valid in a neighborhood of e = 0. The 



result ( |1 . 1| ) then follows by simply expanding in e under the integral sign, eliminating the 
need to fully evaluate any loop integrals. Consequently, in our calculation we never need 
to use the complicated explicit formulas for the e series expansions of M^\e) and M^\e) 
in terms of polylogarithm functions. The only ambiguity introduced by acting with 
is just a polynomial in ln 2 (— t/s), which we show can be fixed up to an additive numerical 
constant by the known infrared and collinear behavior of the two-loop amplitude. 

In order to provide some motivation and context for our calculation we begin in 
section 2 with a more technical review of the difficulties involved in checking (|1 . 1|) and an 
outline of how we propose to overcome them. We present our calculation of the ABDK 
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relation for n = 4 in section 3. Section 4 contains a general argument which proves that 
at any loop order L it is always possible to construct a differential operator which 
generalizes the nice functionality of C^ 2 \ We write down a particularly simple guess for a 
possible choice of based on some very plausible assumptions about the structure of 
higher loop amplitudes. 



2. Advanced Introduction and Motivation 

The ABDK conjecture is highly nontrivial because infrared divergences in 
which start at 0(1/ e 2 ) , imply that it is necessary to know (e) through 0(e 2 ) in order to 
check (T7T) through 0(e ). The proof of (ITT) for n = 4 given in Q proceeds by evaluating 



both sides explicitly in terms of polylogs and finding that they agree up to terms of 0(e). 

Even for n = 4, the direct evaluation of the ABDK relation is difficult for two reasons. 
First of all, both M^\e) and M^\e) are very complicated functions of the kinematic ratio 
x = —t/s. Even using state-of-the-art techniques, such as Mellin- Barnes representations, 
the computation of the e expansions of these two amplitudes requires the resummation of 
an infinite number of poles that give rise to various poly logarithms. To the order required 
to verify these amplitudes involve polylogarithms of degree up to and including 4. 

Secondly, even after M^\e) and M^\e) have both been evaluated to the required order, 
it is necessary to use nontrivial polylog identities to check ( |1 . 1| ) explicitly order by order. 

The goal of this paper is to provide a far simpler method for studying iterative relations 
which we hope may shed some light on the general structure of these relations and may 
be useful to test the conjectures in cases where their status is not currently known. Our 
direct proof relies on three basic ingredients. 

First of all, we use the fact that loop amplitudes can be expressed in terms of Mellin- 
Barnes integral representations (a detailed introductory treatment can be found in the 
book ||22|| ). These have the nice feature that the x dependence can be isolated in a simple 
manner. In particular, any four-particle loop amplitude can be written as an integral of 
the form J^°° dy (—x) y F(e, y) for some function F(e, y). Often it is not very hard to find 
the function F(e,y), rather it is the final integral over y which is exceedingly difficult to 
evaluate. It is clearly tempting to try to collect all of the terms in ( p..l| ) under a single y 
integral and check the ABDK conjecture by expanding in e 'under the integral sign.' In 



other words, we might like to look at the inverse Mellin transformB of (|1 . 1|) with respect 
to the variable — x. 

However, this is usually not possible. The reason is that Mellin-Barnes representations 
of the form J^°° dy (—x) y F(e, y) are generally only valid for e in an open set which does 
not contain e = 0. In order to make sense of statements such as (|1 . 1| ) , which involve a 
series expansion around e = 0, it is necessary to analytically continue the amplitudes in e. 
Unfortunately, as one analytically continues in e towards e = one frequently finds that 
along the way the y integration contour will hit poles in F(e,y). In crossing these poles 
one picks up residue terms which no longer have y integrals. In this paper we will refer to 
any such term which spoils our ability to collect everything under a single y integral and 
then expand around e = as an obstruction. 

We overcome this difficulty with a second ingredient: linear differential operators 
which can be used to eliminate all obstructions, or at least push them off so that they 
are 0(e) and can be ignored. At two loops, the simplest example of an operator which 
accomplishes this goal is (x-^) 5 . However we prefer to use the more elegant choice 

C<$ = (x4--ey (x4-+e)* (2-1) 



dx J \ dx 

which has several attractive features, even though it is of higher degree!. Most importantly, 
we will see that kills all obstructions exactly, to all orders in e. In fact we will see 
that the choice of is completely natural from the Mellin-Barnes integral representation 
of the amplitudes: acting with explicitly removes a number of poles in F(e,y). In 
contrast, acting with (x4-) b does not kill any of the poles; the residues of all of the original 
poles must still be calculated and shown to be O(e) before one can be sure that they can 
be ignored. Also as a minor bonus we note that conveniently preserves the manifest 
x 1/x symmetry of four-particle amplitudes. 

Armed with jC^ we present an elementary proof of the identity 

Mf (6)-i(M«(^ 2 



£ (2) 



0(e), (2.2) 



1 Actually, since the contour for y is taken from — ioo to +ioo, it is more appropriately viewed 
as a Fourier transform with respect to ln(— x) rather than a Mellin transform with respect to —x. 

2 In fact there is even a degree 4 operator (x-^ — e) 2 (x-^ + e) 2 which kills many, but not all, 
of the obstructions, leaving a remainder which is just an innocent number plus 0(e). This operator 
would be sufficient for verifying the ABDK conjecture, although it complicates the analysis slightly 
because one has to keep track of the residues which are not completely killed but only injured. 
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without needing to evaluate the y integral. The proof involves only the evaluation of a finite 
(and small!) number of residues and the use of a single Barnes lemma integral!!. Clearly 
( |2.2|) is a weaker statement than because we can add to M4 (e) — h(M^' (e)) 2 any 
function K(e,x) in the kernel of C^ 2 \ It is clear upon inspection that this kernel consists 
of just a simple function of ln 2 (— x) — therefore, all polylogarithms in the ABDK relation 
are unambiguously fixed already by ( |2.2| ). This involves nontrivial cancellation between 
polylogarithms appearing in the two terms, but we promise the reader that not a single 
polylogarithm will appear in this paper. For us all of the cancellation in ( [2.2| ) occurs under 
the y integral. 

The final ingredient in our proof is the known infrared and collinear behavior of the 

(2) 

two-loop amplitude M\ . We will show that these constraints fix the ambiguity discussed 
above to be of the form K = f^\t)M^ ] (2t) +C + 0(e), where C is a numerical constant 
which our methods do not determine. In this particular case we happen to know from the 
work of 0] that C = — 7r 4 /72, but in general it would probably not be terribly difficult to 
determine C — once it is known to be independent of x — by performing the Mellin-Barnes 
integrals numerically at any convenient value of x (see for example [p3[,|2~4}| ) . 



3. The Two-Loop Calculation 

In this section we begin by proving the identity 



£(2) 



0(e), 



(3.1) 



(2) 

and then use the known infrared and collinear behavior of M\ to fix the ambiguity arising 
from the differential operator C^ 2 \ Of course, one could verify this identity by directly 
substituting the known expressions for Mf } (e) and \{M^\e)f through 0(e) and hitting 
them with C^ 2 \ Our goal in this section is to demonstrate that it is straightforward to 
prove ( |3.1| ) directly at the level of the Mellin-Barnes integrand, without having to first fully 
evaluate the Mellin-Barnes integral in terms of polylogarithms. Along the way we will see 
the special feature of the particular choice of the differential operator which makes 
this possible. 



3 Barnes lemmas are integrals of products of T functions that give rise to more T functions 
and their derivatives. 
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3.1. The |(M 4 (1 } ) 2 Term 

Let us first look at the second term in ( |3.1| ). Taking a convenient expression for 
the one-loop box integral from (7) of ]25|] (see also P^ , P7| ), we find the following simple 
Mellin-Barnes representation for the one-loop amplitude: 



2 (sty/ 2 T(-2e) 



I ¥~ (-x) z+£/2 T(l + e + z)T 2 (z)T 2 (-e - z)T{l - z). 
J 2m 

(3.2) 

The contour for the z integral runs from —zoo to +ioo and passes to the right of all poles 
of the two r(- ■ • + z) functions and to the left of all poles of the two V(- • ■ — z) functions. 
The contour can be taken to be a straight line parallel to the imaginary axis as long as the 
arguments of all V functions have a positive real part@. This is only possible if Re(e) < 0, in 
which case we can take < Re(z) < — Re(e). Note that M^\e) is symmetric under s <-> t, 
or equivalently x 1/x, as can be seen by making the change of variable z — > —z — e. 
Using (|3.2| ) we can write 



w — z 



(3.3) 



1 (,Ai) f nA 2 _ 1 e 2e7 f dz dw 

2 I 4 {e) ) ~ 8(styr*(-2e) J 2m2m [ ~ X ' 

x T(l + e + z)T 2 (z)T 2 (-e - z)T(l 
x T(l + e + w)T 2 (w)T 2 (-e - w)T(l - w). 

Now make the change of variables 

z = u-^y, w = u + ^y (3.4) 

to consolidate the x dependence into the factor (— x) y . This is convenient because acting 
with the differential operator simply inserts a factor of (y 2 — e 2 ) 3 into the integral. 
Note that we hold the product st fixed while differentiating with respect to x = —t/s. 
Therefore 

1 e 2e7 f du dy 



1 /^di, x\ 2 1 e 7 f du dy , . 
- M| 1} 6 = - / —-^F(u,y) 3.5 

2 V / 8 {st) e L z (— 2e) J 2m 2m 



4 We remind the reader that this requirement comes about in the following way. Many of 
the F functions appearing in Mellin-Barnes integral representations originate from integrals of 
the form f da a a_1 (l — a) 6-1 = T(a)T(b)/r(a + b) over Feynman parameters a. This integral 
only converges for Re(a),Re(6) > 0. Other T functions come directly from the Mellin-Barnes 
representation of quantities such as (X + Y)~ u ; the arguments of the T functions arising in this 
way must also have arguments with positive real parts. 
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where 

F{u, y) = (-x)y(y 2 - e 2 ) 3 r(l + e + u - \y)Y{u - \y) 2 Y{-t -u+ §y) 2 r(l -u + \y) 
xT(l + e + « + \y)T{u + \y) 2 Y{-e -u- \y) 2 Y(l -u-\y). 

(3.6) 

Let us now reviewi the procedure for manipulating Mellin-Barnes integrals such as 
( p.5[ ). The formula (|3.5|) defines a meromorphic function of e for e in an open set which 
does not contain e = 0. As we take e — > 0, the integration contours become pinched 
until at e = they are forced to run through some poles of the Y functions. In order to 
construct a formula for (AQ (e)) 2 which is valid in an open neighborhood of e = we 
must analytically continue ( |3.5|) in e. 

For definiteness, let us choose Re (it — hy) > Re(u + hy) > 0. The contour C for the u 
and y variables must be such that the argument of each Y function has a positive real part. 
This can be satisfied as long asi Re(e) < Re(— u — hy) < 0. With this choice, the first 
pole that hits one of the integration contours as we try to take e — > is the first pole of 
r(— e — u + hy)- Passing the pole through this contour produces two terms: the first term 
is the same integral as in ( |3.5| ), but with a contour C that now passes to the right of the 
first pole of r(— e — u + hy), and the second term is the residue of F(u, y) at u = — e + ~y. 
The residue comes with a minus sign because the contour passes across the pole from left 
to right, leaving a clockwise integral around the pole. Therefore we have 

/ = I ~ I ^->- + i»nu,y). (3.7) 

Jq ZTll LTXl J c , Z7TZ Z7TZ Jq Illl 2 " 

Now when Re(e) becomes larger than Re(— u + hy), the contour C can be taken to be the 
straight line C. Therefore we replace C by C in (|3.7|) to write 



/ ^^-F{u,y)^ [ ^-.^-.F{u,y) - [ ^Res u= _ eHy F(u,y). (3.8) 
J c 2m 2m J c 2m 2m J c 2m 2y 

The notation is a reminder of the logic here: it would be incorrect to write =, since that 
would require the second term on the right-hand side to be zero. The left- and right-hand 



5 See [26,28] for pioneering work on how this is implemented at two loops. A very convenient 
program which automates these kinds of manipulations has recently been made available by 



M. Czakon f2J. 

6 There is also a lower bound on Re(e). We do not worry about this since we are always 
interested in pushing e towards 0, not away from it. We can imagine starting out with Re(e) 
infinitesimally less than Re(— u — 
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sides of (3^8) represent the same meromorphic function of e. The left-hand side is valid 
in a neighborhood of Re(e) < Re(— u + ^y) while the right-hand side defines the analytic 
continuation of the left-hand side for a neighborhood Re(e) > Re(— u + \y). 

The only remaining pole which we encounter as we take e to zero is the first pole of 
-e — \y. This gives 

du dy 

F(u,y) 

(3.9) 



T(— e — u — at u 
du dy 



c 27ri 2-rci 



F(u,y) 



c 2ixi 2-rci 
dy 



2txi 



-Res, 



iv F ( u >y) - 



dy 
2m 



Res, 



F{u,y), 



with the left-hand side defined for Re(e) < Re(— u + \y) and the right-hand side defined 
for Re(e) > Re(— u — \y)- There are no more contours standing in the way of taking 
e — > 0, so the three terms on the right-hand side of (|3.9| ) can therefore be evaluated in an 
open neighborhood of e = by simply making a power series expansion in e under the 
y integral. The first term in ( |3.9|) (the double integral term) is manifestly 0(1) because 
F(u,y) itself is. When we remember the prefactor in front of the integral ( |3.5| ), which 
contains an explicit l/T(—2e), we see that this term only contributes to ( |3.5| ) at 0(e), so 
we can ignore it. All we have to do is evaluate the two residues on the second line of 
and expand the result (including the prefactor in (|3.5|) ) through O(e ), which gives 



M«(e) 



dy , 

(— X 

2ni K 

2 

€ 



Vr(i-y)r(-y) 2 r(y)* 2 r(i + y ) 



-- e +y(f (2) (y)-f (2) (-y) 



(3.10) 



+ 0(e), 



where is the function defined in ( |1.2| ). It is very intriguing to see this function, which 
is related to the infrared and collinear behavior of two- loop amplitudes, emerge in an 
interesting way from a one-loop amplitude squared. We follow the notation of [22| in using 



r* as a reminder that the contour for the y integral passes the first pole of the indicated V 
function T 2 (y) on the wrong side. The appropriate contour for ( |3.10| ) is — 1 < Re(y) < 0, 
reflecting a choice we made in setting up the original u and y contours in ( |3 . 5|) . 

It may appear that the application of the differential operator C^ 2 \ and the corre- 
sponding insertion of the factor (y 2 — e 2 ) 3 into the integral, was not important in this 
analysis. In fact it plays a crucial role in removing a third-order obstruction. Without 
including this factor in F we would have found 

2(-x)T(-e) 4 r(l + e) 4 



-Res u= _ e+ i y F(u,y) 



+ o((e- y y 



(3.11) 



which contains an obstruction to taking e — > in the form of a triple pole at y = e. A 
factor of (y — e) 3 is needed to kill this obstruction. The necessity of the other factor (y + e) 3 
only becomes apparent in the analysis of M\ (e), to which we now turn our attention. 



3.2. TheM {2) Term 



In the previous subsection we reviewed in detail the procedure for manipulating Mellin- 
Barnes integrals. In this section we will calculate C^M^^e) by the same procedure, 
although we will not show each step in as much detail. The full two-loop amplitude is 



given by the sum of the two terms [13] 



Mf(e) = - A s 2 tlf\s,t) + ±st 2 li 2 \t,s), (3.12) 

where I\ (s,t) is the two-loop massless scalar box function. We start with a convenient 
Mellin-Barnes representation for the two-loop box function given by (11) of rf: 

1 -2 +r (2)/ A _ 1 e 2e7 f da dz dz x dz 2 , u_ ff+e 



-sHr >(s t) = - - / ^Ljnmm (- x \ 

4 4 v ' ; 2 (st) e r(-2e) J 2ni 2ni 2ni 2ni v ; 

r(-zi - z 2 )T(-z 1 + z 2 )T(2 + e + z + gi) T(q)r(l + e - z - - a) 

r(l - zi - z 2 )T(l - z x + z 2 ) r(l-2e + 2 + ^i) ( 3 - 13 ) 

x r(-l -e-z- z 2 )T(-l -e-z + z 2 )Y{l + z- Zl )T(l + z + z x ) 

x T(-e + z + z 2 + a)T(-z + z x - a)T(-e + z - z 2 + a)T(l - a). 

The full two-loop amplitude ( |3.12|) is obtained by adding together this quantity and the 
same expression with x replaced by 1/x. Our goal is to reduce this to a single integral over 
the variable a, since that is the only variable which sets the x dependence through the 
factor (— x) 1 ~ a+£ . The contour is fixed by requiring that the argument of each function 
has a positive real part. The final result does not depend on the particular choice, but 
intermediate steps can look a little different. For definiteness, we imagine taking the 
contour Re(z) = -2/3, Re(zi) = -1/8, Re(^ 2 ) = -1/16 and Re(tr) = 1/2. 

A simple analysis of ( |3.13| ) shows that all of the poles we hit as we analytically continue 
e to sit at only four different possible values of z 2 , 

z 2 = l + e + z, z 2 = —l — e — z, z 2 =e-a-z, z 2 = -e + a + z. (3.14) 



7 Our formula follows from theirs after the change of variables a = z\ + Z2, ft = Z\ — Z2, 
t = z — Z\. 
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Our first step is therefore to deform the contour past these poles. This leaves us with a 
sum of four residue terms, each of which is only a triple integral, plus the original four- 
fold integral ( |3.13|) , which can now be evaluated in a neighborhood of e = 0. In fact the 
four-fold integral is manifestly O(e) due to the explicit l/r(— 2e) factor sitting in front of 
( |3.13|) . Therefore we only need to follow the four residues ( |3.14j) . Let's look at the first 
residue term, 

-sHl\ 2) ( 8 .t) = -—- ; r / (-X y- a+e 

4 4 v ' ' 2 (styr(-2e) J 2m 2m 2ni K ! 

r(-l - e - z - z x )Y{l + e + z- z 1 )T(2 + e + z + z x ) r(g)r(l + e-z-z 1 -a) 
X r(-e - z - z 1 )T{2 + e + z - zi) r(l - 2e + z + Zl ) 

x r(-2 - 2e - 2z)T(l + z- zi)T(l + z + z x ) 

x T(l + 2z + a)Y(-z + z x - a)r(-l - 2e + a)Y(l - a) 

+ three more residues + 0(e). 

(3.15) 

This expression is now valid for —11/24 < Re(e) < —1/3 (given our specific choice for the 
contours). As we continue taking e closer to zero we encounter more poles. The ones we 
are concerned about are poles which force us to take a residue in a. For example, when 
Re(e) reaches —1/4 we hit a pole at a = 1 + 2e. The residue at this pole is a double integral 
J dz\dz. The presence of this term spoils our goal of trying to perform all manipulations 
under the a integral. 

The dangerous pole at a = 1 + 2e comes from the T(-l - 2e + a) factor in (g7T5|) . 
We could get rid of this pole by inserting a factor of —1 — 2e + a into the integral. But 
this is precisely the factor we would get if we acted on ( |3.15| ) with the differential operator 
—x4^ — e! Applying this differential operator kills the pole completely; the corresponding 
residue is exactly zero, not just zero to 0(e). 

It is a straightforward exercise to continue following the Mellin-Barnes procedure all 
the way to e = 0. It turns out to be necessary to cube the factor (—1 — 2e + a) inside 
the integrand as this kills other obstructions which appear at later steps. This analysis 
motivates consideration of the differential operator (—x-^ — e) 3 . However, we have only 
looked at the term js 2 tl^\s, t). The complete two-loop amplitude (|3.12| ) is the sum of 
this and ^st 2 I^ (£, s), which is the same up to x — > \ jx. The same analysis applied to 

(2) 

this other term suggests that we should also hit M\ (e) with the differential operator 
(x-4- — e) 3 . Taken together, these observations motivate the choice of in ( |2.1|) . 
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Let us now display the final result for acting on the Mellin-Barnes integral ( |3.15|) . 
Each of the triple integral terms in ( [3.15 ) eventually branches into the same triple integral, 
now valid in a neighborhood of e = 0, plus a sum of residue terms which involve double 
or single integrals. All four triple integral terms, and several double integral terms which 
appear later on, turn out to be 0(e), again because of the explicit factor 1/T(— 2e) in 
( |3.15|) . However, this explicit factor 1/F(— 2e) certainly does not mean that everything 
is 0(e) since taking residues can produce explicit singularities which cancel this factor. 
Ultimately we find that only 6 residues contribute through 0(e), 



~sHl?\ Sl t) 



1 



3 2e 7 



2 (styT(-2e) 



da 

— Res 2l=e Res 2 

ZTVl 



-1-e 



Res 



Z 2 = l + £ + Z 



+ 



da 



i-.+ i^Res, 



+ 
+ 



Res^ _ i i rr Res 2 -_i 

da 

Res 7l - 

2ni 



da 
2tH 
da 
2tH 

da 

2tH 



Zl Res 22= i+ e + 
i_ e _ 2l Res 



3 Z2 = 



1-e-z 



(3.16) 



R6s 2l=e Res 2=e _ -Res 2;2=e _ 2 _ ( j 



dz\ 
2tH 

dz\ 
2tvl 



Res 



z = e + Zi — cr 



Res 



Z2=£ — z—a 



Res 



z = e + Zi — cr 



Res 



Z2 = 



-e+a+z 



G + 0(e), 



where G is the quantity appearing inside the integral in (|3.13| ). The signs are easily fixed 
by keeping in mind that a residue of the form Res„ = ...± e comes with a sign of ±1. 

A simple calculation reveals that the second and third line in (|3.16|) are equal to each 
other, as are the fifth and sixth lines. Let us look at these last two terms first since they 
still involve an extra z\ integral. Expanding to 0(e) we find that they contribute 



1 



[sty 



da 
2tH 



-x 



.i-ff 



(i-a) t) r( C T)r(cx-i)r* 2 (i-a 



[p : T(2-a + 2z 1 )T*(2z 1 )r(-2z 1 )T*(a - 1 - 2z x ) + 0(e), 
J 2tci 



(3.17) 



where as before we keep track of the integration contour by using a * to denote a V 
function whose first pole is crossed in the 'wrong' direction. The remaining z\ integral can 
be performed with the help of the generalized Mellin-Barnes identity 

dz 



2ni 



r(Ai + ^)r*(A 2 + z)T(-X 2 - z)T*(X 3 - z) 



T(A 2 + A 3 ) (r(Ai - A 2 ) [^(Ai - A 2 ) - ^(Ai + A 3 )] - T(-A 2 - A 3 )r(Ai + A 3 )), 



(3.18) 
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which follows straightforwardly from (D.2) of ||22|| . 

Combining (|3.17| ) with the first four terms in Q3.16|) , expanded through 0(e ), gives 



£(2) 



— +47rcot(7r(l - a)) - ln(-x) +^(2 - a) - ^(1) 



+ 0{e). 
(3.19) 

To get the full two-loop amplitude ( |3.12| ) we should add ( |3.19| ) to the same quantity 
with x replaced by 1/x. This looks equivalent to leaving (— x) 1_<J alone and replacing 
a — ► a' = 2 — a (and ln(— x) — > — ln(— x)) inside the integral, but there is a subtlety. The 
contour for a was chosen to run along Re(cr) = -|, which would place the contour for a' 
at |. We would like to rename a' back to a and combine both terms under a single a 
integral. In order to do this we must check that we don't cross any poles in taking a' from 
| — > Fortunately the high power (1 — cr) 6 is more than adequate to kill the pole at 
a = 1, so there is no problem. 

Finally we conclude that 



£ (2) 



Mf(e) 



1 



(sty 



^ {-xf-'il - a) 6 I»r 2 (a - l)r* 2 (l - a)T(2 - a) 
2m 



+ V(2-t7)-2V>(l)+V(<r) 



(3.20) 



+ 0(e). 



The change of variable a — > 1 — y and the identity 



t/>(1 -y) + ^(1 + y)= ^(y) + ^(-y) 



(3.21) 



precisely transform ( 3.20 ) into (|3 .10 ). This concludes the elementary proof that 



£ (2) 



0(e). 



(3.22) 



Note that while M"f } (e) and ±(Mf } (e)) 2 are both individually C(e" 4 ), ( CT ) and (Q) 
indicate that they become C(e _1 ) after being hit with the differential operator £( 2 **. This 
fact is evident when one looks at the explicit formulas 



( S t) e Mf( e ) = ^-l 



Iln 2 (-s) + ^ 2 



{stf\ (Mi 1} (e) 



2 _ 2 

~~ 74 



1 



1 



ln 2 (-x) + -tt' 



+ 0(e- 1 ] 



(3.23) 
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There is at least one slightly intriguing aspect of the calculation in this subsection. In 
the calculation of |4|] it was remarked that there is a non-trivial cancellation of terms, in- 

(2) (2) 

volving the use of polylogarithm identities, between the two integrals I\ (s, t) and I\ (t, s) 
which contribute to the two-loop amplitude. In our calculation this cancellation manifests 
itself under the a integral as 

cot(7r(l - a)) + cot(yr(a - 1)) = 0. (3.24) 

3.3. Fixing The Ambiguity 

The result QTgg) implies that 

Mf\e) - \ (M| 1} (e)) 2 = ^K(x,e) + 0(e), (3.25) 

where K (x, e) is a conveniently normalized undetermined function in the kernel of C^ 2 K To 
conclude the check of the ABDK relation (|1 . 1|) we have to study this kernel. Now K(x, e) 
has to be invariant under x — > \ jx (since the left-hand side of (|3.25 ) is), and it is easy to 



see that the most general function annihilated by which respects this symmetry is 

K(x,e)=A(e) ((-x) e + (-x)~ £ ) 

+ B(e)ln(-x) ((-x) e - (-x)" e ) (3.26) 
+ C(e) \n\-x) {{-xf + (-x)" e ) , 

where A(e), B(e) and C(e) are arbitrary functions of e but independent of x. Note that 
( ^.26| ) is a function of ln 2 (— x) only. 

By inspecting the original relation ( |1 . 1| ) it is clear that an important consistency check 
of our method is to find that (st) e M4 (2e) can be expressed in the form ( |3.26| ) up to terms 
of order e. Indeed, a little algebra shows that 



( S t) e M| 1} (2e) 



4e 2 + 3 



{(- x y + (-x)- e ) + !^L^ ((- x y - (-x)- e ) + 0( e ). 



(3.27) 

Equivalently, it is of course also possible to check directly that M^\2e) = O(e). Since 
we have verified that the quantity (si) e iif| 1 ' ) (2e) lies in the kernel of (to 0(e)), we can 
freely absorb this into the ambiguity in K(x, e) to rewrite ( |3.25|) in the equivalent form 

Mf\e) - ~ (Mi 1} (6)) 2 - f^M { l\2e) = ^K(x, e) + 0(e). (3.28) 
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Let us now take a look at what we can say about the remaining ambiguity K(x : e) 
from the infrared singularity structure of two-loop amplitudes. This has been studied in 
QCD and can be applied to M = 4 SYM as well. The known behavior implies that the 
quantity appearing on the left-hand side of ( ft.28|) is free of IR divergences. This implies 
that K (x, e) must not have any poles in e. Looking at the most general form ( |3.26 ) of 
K(x, e), it is easy to see that the only function compatible with the infrared behavior is 

' K(x,e) = C + E\n 2 (-x) +Fln 4 (-x) +0{e), (3.29) 



(st) 

where C, E and F are numerical constants. 

In fact it is easy to reduce the ambiguity slightly (|3.29| ) by using a lower-order differ- 
ential operator instead of £( 2 >. For example, the operator (x-^ + e) 2 {x-^ — e) 2 , which we 
mentioned in section 2, kills most, but not all, of the obstructions, leaving as remainders 
just numerical constants (which can be absorbed into C in ( |3.29| )). One could even con- 
sider just the operator (x^-) 4 , which does not completely kill any obstruction, but reduces 
all of them to numbers plus 0(e). The use of these operators would have complicated the 
analysis of section 2, since we would have to keep track of all of these unkilled obstructions. 
Their only advantage is that these operators have smaller kernels than C^ 2 \ so one can 
use them to argue that F in ( |3.29| ) must be zero. 

An independent argument which reduces the ambiguity ( |3.29| ) even further involves 

(2) 

considering the collinear limit of the two-loop amplitude (e). The known behavior 
P,p|-p!2| implies that the quantity on the left-hand side of ( |3.28|) must be finite as x — > 
and x — > oo, which immediately fixes E = F = 0. In conclusion, we have verified the form 
of the ABDK relation 

M {2) (e) = \ (m| 1} (e)) 2 + f {2) (e)Mi 1} (2e) + C + 0(e) (3.30) 

up to an overall additive numerical constant. 



4. Higher Loops 

In this section we would like to tie together the threads which have been weaving 
around throughout the analysis of sections 2 and 3 into a coherent picture. We show 
that the methods used in this paper can be generalized for n = 4 particle amplitudes at 
L > 2 loops in a completely straightforward manner. We first prove that it is always 
possible to find a generalization of the differential operator which naturally removes 
all obstructions to combining all of the terms in the L-loop four-particle iteration relation 
inside a single integral. In subsection 4.2 we will write down a particularly simple form for 
£( L ) which might be a natural guess to try. 
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4-1. General Constructions 

The multiloop iteration relations proposed in [BJ take the form 



Mk L \e) =xM[Mg)(e)] + f {L \e)M^(Le) +C™ +0(e), (4.1) 

where are numerical constants, are functions of e only, and is a polynomial 
in the lower loop amplitudes which is conveniently summarized in the expression 



Xf)[M«]=Mf)-ln l + 5>^4° 



2 = 1 



(4.2) 

a L term 



For L = 2 we recover ( |1 . 1| ) while for L = 3 we have for example 

M< 3 >(e) = ~ (M^ie)) 3 + M^\e)M^\e) + /( 3) (e)M«(3e) + + O(e). (4.3) 

This relation has been shown to be correct for n = 4 by explicit evaluation of the partici- 
pating amplitudes 0. 

We would like to collect all of the terms in ( |4.3|) or its higher L generalizations, for 
n = 4 particles, under a single Mellin integral of the form 

dy (-x) y F(e,y) (4.4) 

and then to prove the iteration relations by expanding in e under the integral sign. We 
found that in general this was not possible due to obstructions which prevent us from 
taking e — > inside the integral. The most general possible obstruction is a residue of the 
form 

Res v=s(c) [(-x)"F(e,y)] (4.5) 

where g(e) is determined by the precise arguments of the V functions appearing in F(e, y). 
In general it will be a linear function of e. (For example, the obstruction we discussed 
below ( 3.15|) corresponds to g(e) = — e when the a variable is transcribed to the y variable 



used in this section.) The crucial feature of the expression ( |4.5|) is that while the resulting 
e dependence can be very complicated, the x dependence is simple. Let us suppose that 



the pole where we are taking the residue is a pole of order k. Then the residue ( |4.5| ) will 
evaluate to 

(-x)*( e >P fc (ln(-x),e), (4.6) 
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where is some polynomial in In (—a:) of degree k. The coefficients of that polynomial 
might have very complicated e dependence. A differential operator which exactly kills the 
residue ( }4.6| ) is 



x-- 9(e) j . (4.7) 

In (|4.5|) we analyzed the contribution from a single obstruction. In a general amplitude 
(or product of amplitudes, such as appear in the iteration relations) there will be many 
such obstructions. Obviously it is always possible to construct a differential operator which 
removes all of these obstructions by taking the product of all of the individual differential 
operators of the form (|4.7|) — one for each obstruction. 



Now suppose that we would like to test a conjectured iteration relation for M± (e) 
using some differential operator constructed in the manner we have just described. 
By construction, applying to the iteration relation will allow us to combine all of 
the terms under a single y integral, and then to verify the iteration relation by expanding 
the integrand through 0(e ) inside the integral. The ambiguity introduced by acting 
with consists precisely of terms of the form ( |4.6| ), which span the kernel of (by 
construction). Since g{e) is always a linear function of e, we can expand fl4.6| ) in e to see 
that the only things which can appear are powers of ln(— x) (actually, the ambiguity must 
be a polynomial in ln 2 (— x), due to the x — > 1/x symmetry). There can never be anything 
as complicated as a polylogarithm function — those are uniquely fixed even after is 
applied. 

In summary, we have demonstrated in this subsection that it is always possible to 
construct a differential operator which, when applied to a conjectured iterative re- 
lation for M4 (e), allows the relation to be checked without explicitly evaluating the y 
integral. Moreover all of the complicated generalized polylogarithm functions are com- 
pletely fixed by this procedure; the only ambiguity introduced by is a polynomial in 
ln 2 (— x). It might be possible to argue away some of this ambiguity by a more careful 
general argument. 

4-2. Specific Conjectures 

An example of an L-loop generalization of that would be natural to try is 

mv\ f d L \ L+1 ( d L \ L+1 , A . 

£( L ) = x — + -e) x- e 4.8) 

V dx 2 J V dx 2 J K 1 
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The form of this guess is based on some observations of the particular structure of the 
Mellin-Barnes representations that one encounters at one, two, and three loops. 

The first observation is that obstructions tend to be of the form (|4.6|) where in fact 
g(e) = ±Le/2. We are not aware of any case where this is not true, so we conjecture that 
this might be a general feature. If true, it would mean that the operator ( |4.8|) kills all 
obstructions exactly, to all orders in e. Even if (|4.8 ) does not kill all obstructions exactly, 
it is clear that it is of sufficiently high order in j£ to push all obstructions to 0(e 2 ) or 
higher, where we do not care about them. 



5. Summary and Outlook 

In this paper we have proposed to use certain simple linear differential operators as 
effective tools for studying the iterative structure of planar four-particle L-loop MHV am- 
plitudes in J\f = 4 SYM. A key ingredient in our analysis is played by Mellin-Barnes integral 
representations of loop amplitudes, and by what we call 'obstructions.' Obstructions are 
terms which appear as we analytically continue Mellin-Barnes integrals towards the region 
around e = 0; their presence implies that it is impossible to collect everything under a sin- 
gle integral which simultaneously preserves the x-dependence in the simple factor (— x) y 
and also admits a series expansion in e under the y integral. We showed that obstructions 
are always functions of ln 2 (— x) and may be killed by the application of simple differential 
operators. The differential operator may be chosen to be simply x£ to a sufficiently high 
power if one is only interested in pushing the appearance of obstructions to a sufficiently 
high order in e, or a slightly more complicated operator if one wishes to kill them exactly. 
The operator defined in ( j2.1| ) was proven to kill all obstructions appearing in the 
two-loop ABDK iterative relation (|1.1|) . We argued that the generalization guessed 
in Q4.8| ) might achieve this task at L loops. 

The advantage of removing obstructions is obvious: by getting rid of them, we can 
directly study the inverse Mellin transform of an iterative relation in an expansion around 
e = 0. In other words we have no need to evaluate the final, and often exceedingly 
complicated, y integral explicitly in terms of generalized polylogarithms. The resulting 
simplicity is nicely exhibited by our elementary verification of the structure of the two- 
loop ABDK relation ( |1 . 1|) for n = 4. The proof of ( |3.22|) for required nothing more than 



evaluating a small number of residues and using the single Mellin-Barnes integral (|3.18| ). 
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Fortunately, the price we pay for this beautiful simplification is not too high: by 
construction, the only ambiguities introduced by the differential operators we study are 
polynomials in ln 2 (— x). The complicated generalized polylogarithms are unambiguously 
fixed. It is likely that the remaining ln 2 (— x) ambiguity can in general be fixed by the 
known infrared and collinear behavior of L-loop amplitudes. 

It is natural to suppose that these techniques can also be generalized to n > 4 particle 
amplitudes. In general these admit Mellin-Barnes representations where we isolate the 
dependence on the various independent kinematic variables, 

/ ^(h-k 2 r [ ^(h-k 3 )y*---F(e ]yi ,y 2 ,---). (5.1) 
J 2m J 2m 

There will now be a more complicated picture involving complete obstructions and partial 
obstructions in some subset of the variables, but it should nevertheless be possible to 
construct partial differential operators in several variables which might prove useful in 
studying the iterative structure of these amplitudes. 

Another interesting direction might be to use the converse mapping theorem to com- 



pute the asymptotic behavior of Feynman integrals [|30| as t/s — ► and as t/s — > oo in 
order to study possible iterative relations in cases where proving the full structure is not 
within reach!. Particularly interesting would be to study the n = 4 four-loop case. 
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